Deepening genomic sequences of 1081 Gossypium hirsutum accessions reveals novel SNPs and haplotypes relevant for practical breeding utility.

Fiber quality is a major breeding goal in cotton, but phenotypically direct selection is often hindered. In this study, we identified fiber quality and yield related loci using GWAS based on 2.97 million SNPs obtained from 10.65× resequencing data of 1081 accessions. The results showed that 585 novel fiber loci, including two novel stable SNP peaks associated with fiber length on chromosomes At12 and Dt05 and one novel genome regions linked with fiber strength on chromosome Dt12 were identified. Furthermore, by means of gene expression analysis, GhM_A12G0090, GhM_D05G1692, GhM_D12G3135 were identified and GhM_D11G2208 function was identified in Arabidopsis. Additionally, 14 consistent and stable superior haplotypes were identified, and 25 accessions were detected as possessing these 14 superior haplotype in breeding. This study provides fundamental insight relevant to identification of genes associated with fiber quality and yield will enhance future efforts toward improvement of upland cotton.

Identification of Crucial Modules and Genes Associated with Bt Gene Expression in Cotton.

The expression of Bacillus thuringiensis (Bt) toxins in transgenic cotton confers resistance to insect pests. However, it has been demonstrated that its effectiveness varies among cotton cultivars and different tissues. In this study, we evaluated the expression of Bt protein in 28 cotton cultivars and selected 7 cultivars that differed in Bt protein expression for transcriptome analysis. Based on their Bt protein expression levels, the selected cultivars were categorized into three groups: H (high Bt protein expression), M (moderate expression), and L (low expression). In total, 342, 318, and 965 differentially expressed genes were detected in the H vs. L, M vs. L, and H vs. M comparison groups, respectively. And three modules significantly associated with Bt protein expression were identified by weighted gene co-expression network analysis. Three hub genes were selected to verify their relationships with Bt protein expression using virus-induced gene silencing (VIGS). Silencing GhM_D11G1176, encoding an MYC transcription factor, was confirmed to significantly decrease the expression of Bt protein. The present findings contribute to an improved understanding of the mechanisms that influence Bt protein expression in transgenic cotton.

Identification and Expression Analysis of EPSPS and BAR Families in Cotton.

Weeds seriously affect the yield and quality of crops. Because manual weeding is time-consuming and laborious, the use of herbicides becomes an effective way to solve the harm caused by weeds in fields. Both 5-enolpyruvyl shikimate-3-phosphate synthetase (EPSPS) and acetyltransferase genes (bialaphos resistance, BAR) are widely used to improve crop resistance to herbicides. However, cotton, as the most important natural fiber crop, is not tolerant to herbicides in China, and the EPSPS and BAR family genes have not yet been characterized in cotton. Therefore, we explore the genes of these two families to provide candidate genes for the study of herbicide resistance mechanisms. In this study, 8, 8, 4, and 5 EPSPS genes and 6, 6, 5, and 5 BAR genes were identified in allotetraploid Gossypium hirsutum and Gossypium barbadense, diploid Gossypium arboreum and Gossypium raimondii, respectively. Members of the EPSPS and BAR families were classified into three subgroups based on the distribution of phylogenetic trees, conserved motifs, and gene structures. In addition, the promoter sequences of EPSPS and BAR family members included growth and development, stress, and hormone-related cis-elements. Based on the expression analysis, the family members showed tissue-specific expression and differed significantly in response to abiotic stresses. Finally, qRT-PCR analysis revealed that the expression levels of GhEPSPS3, GhEPSPS4, and GhBAR1 were significantly upregulated after exogenous spraying of herbicides. Overall, we characterized the EPSPS and BAR gene families of cotton at the genome-wide level, which will provide a basis for further studying the functions of EPSPS and BAR genes during growth and development and herbicide stress.

Lysine 2-Hydroxyisobutyrylation- and Succinylation-Based Pathways Act Inside Chloroplasts to Modulate Plant Photosynthesis and Immunity.

Crops must efficiently allocate their limited energy resources to survival, growth and reproduction, including balancing growth and defense. Thus, investigating the underlying molecular mechanism of crop under stress is crucial for breeding. Chloroplasts immunity is an important facet involving in plant resistance and growth, however, whether and how crop immunity modulated by chloroplast is influenced by epigenetic regulation remains unclear. Here, the cotton lysine 2-hydroxyisobutyrylation (Khib) and succinylation (Ksuc) modifications are firstly identified and characterized, and discover that the chloroplast proteins are hit most. Both modifications are strongly associated with plant resistance to Verticillium dahliae, reflected by Khib specifically modulating PR and salicylic acid (SA) signal pathway and the identified GhHDA15 and GhSRT1 negatively regulating Verticillium wilt (VW) resistance via removing Khib and Ksuc. Further investigation uncovers that photosystem repair protein GhPSB27 situates in the core hub of both Khib- and Ksuc-modified proteins network. The acylated GhPSB27 regulated by GhHDA15 and GhSRT1 can raise the D1 protein content, further enhancing plant biomass- and seed-yield and disease resistance via increasing photosynthesis and by-products of chloroplast-derived reactive oxygen species (cROS). Therefore, this study reveals a mechanism balancing high disease resistance and high yield through epigenetic regulation of chloroplast protein, providing a novel strategy to crop improvements.

Preparation of chitin/MXene/poly(L-arginine) composite aerogel spheres for specific adsorption of bilirubin.

Currently, hemoperfusion is clinically the most rapid and effective treatment for removing toxins from the blood. The core of hemoperfusion is the sorbent inside the hemoperfusion device. Due to the complex composition of the blood, adsorbents tend to adsorb substances such as proteins in the blood (non-specific adsorption) while adsorbing toxins. Hyperbilirubinemia is caused by excessive levels of bilirubin in the human blood, causing irreversible damage to the patient's brain and nervous system, and even leading to death. High adsorption and high biocompatibility adsorbents with specific bilirubin adsorption are urgently needed to treat hyperbilirubinemia. Herein, poly(L-arginine) (PLA) which can specifically adsorb bilirubin, was introduced into chitin/MXene (Ch/MX) composite aerogel spheres. Ch/MX/PLA prepared by supercritical CO2 technology had higher mechanical properties than Ch/MX and can withstand 50,000 times its own weight. The in vitro simulated hemoperfusion test showed that the adsorption capacity of Ch/MX/PLA was as high as 596.31 mg/g, which was 15.38 % higher than that of Ch/MX. Binary and ternary competitive adsorption tests showed that Ch/MX/PLA also had good adsorption capacity in the presence of a variety of interfering molecules. In addition, hemolysis rate testing and CCK-8 testing confirmed that Ch/MX/PLA had better biocompatibility and hemocompatibility. Ch/MX/PLA can meet the required properties of clinical hemoperfusion sorbents and has the ability to produce mass production. It has good application potential in the clinical treatment of hyperbilirubinemia.

Transcriptome, Ectopic Expression and Genetic Population Analysis Identify Candidate Genes for Fiber Quality Improvement in Cotton.

Comparative transcriptome analysis of fiber tissues between Gossypium barbadense and Gossypium hirsutum could reveal the molecular mechanisms underlying high-quality fiber formation and identify candidate genes for fiber quality improvement. In this study, 759 genes were found to be strongly upregulated at the elongation stage in G. barbadense, which showed four distinct expression patterns (I-IV). Among them, the 346 genes of group IV stood out in terms of the potential to promote fiber elongation, in which we finally identified 42 elongation-related candidate genes by comparative transcriptome analysis between G. barbadense and G. hirsutum. Subsequently, we overexpressed GbAAR3 and GbTWS1, two of the 42 candidate genes, in Arabidopsis plants and validated their roles in promoting cell elongation. At the secondary cell wall (SCW) biosynthesis stage, 2275 genes were upregulated and exhibited five different expression profiles (I-V) in G. barbadense. We highlighted the critical roles of the 647 genes of group IV in SCW biosynthesis and further picked out 48 SCW biosynthesis-related candidate genes by comparative transcriptome analysis. SNP molecular markers were then successfully developed to distinguish the SCW biosynthesis-related candidate genes from their G. hirsutum orthologs, and the genotyping and phenotyping of a BC3F5 population proved their potential in improving fiber strength and micronaire. Our results contribute to the better understanding of the fiber quality differences between G. barbadense and G. hirsutum and provide novel alternative genes for fiber quality improvement.

SEP-like genes of Gossypium hirsutum promote flowering via targeting different loci in a concentration-dependent manner.

SEP genes are famous for their function in the morphological novelty of bisexual flowers. Although the diverse functions of SEP genes were reported, only the regulatory mechanisms underlying floral organ development have been addressed. In this study, we identified SEP-like genes in Gossypium and found that SEP3 genes were duplicated in diploid cotton varieties. GhSEP4.1 and GhSEP4.2 were abundantly transcribed in the shoot apical meristem (SAM), but only GhSEP4.2 was expressed in the leaf vasculature. The expression pattern of GhSEPs in floral organs was conserved with that of homologs in Arabidopsis, except for GhSEP2 that was preponderantly expressed in ovules and fibers. The overexpression and silencing of each single GhSEP gene suggested their distinct role in promoting flowering via direct binding to GhAP1 and GhLFY genomic regions. The curly leaf and floral defects in overexpression lines with a higher expression of GhSEP genes revealed the concentration-dependent target gene regulation of GhSEP proteins. Moreover, GhSEP proteins were able to dimerize and interact with flowering time regulators. Together, our results suggest the dominant role of GhSEP4.2 in leaves to promote flowering via GhAP1-A04, and differently accumulated GhSEP proteins in the SAM alternately participate in forming the dynamic tetramer complexes to target at the different loci of GhAP1 and GhLFY to maintain reproductive growth. The regulatory roles of cotton SEP genes reveal their conserved and diversified functions.

Freeze-drying platforms design for batch fabrication of Haversian system mimicking scaffolds with enhanced osteogenesis.

The Haversian system is one of the most important pathways to repair bone defects, and it is the basic guarantee for the repair of bone defects, which means that the formation of the Haversian system indicates repairing of the defects. The integration of structure and function for tissue engineering scaffolds is of great importance in mimicking native bone tissue. However, in contrast to the increasing demands, how to rapidly prepare various sizes of such Haversian system mimicking scaffolds in batch becomes a major challenge. In this study, we designed three types of platforms with different sizes in combination with the freeze-drying approach. Chitosan/type I collagen composite materials were used to study the structure, morphology, and performance of the production, and the effects of the controlled architecture on osteogenesis. Results showed that the physicochemical effects of the mass fabricated scaffolds of various sizes met the requirements of bone repair materials. In addition, the scaffolds had good cytocompatibility and excellent in vivo bone repair performance, which have potential clinical applications.

Development and Utilization of Functional Kompetitive Allele-Specific PCR Markers for Key Genes Underpinning Fiber Length and Strength in Gossypium hirsutum L.

Fiber length (FL) and fiber strength (FS) are the important indicators of fiber quality in cotton. Longer and stronger fibers are preferred for manufacturing finer yarns in the textile industry. Functional markers (FMs) designed from polymorphic sites within gene sequences attributing to phenotypic variation are highly efficient when used for marker-assisted selection (MAS) in breeding superior varieties with longer FL and higher FS. The aims of this study were to develop FMs via kompetitive allele-specific PCR (KASP) assays and to validate the efficacy of the FMs for allele discrimination and the potential value in practice application. We used four single-nucleotide polymorphism markers and 360 cotton accessions and found that two FMs, namely, D11_24030087 and A07_72204443, could effectively differentiate accessions of different genotypes with higher consistency to phenotype. The appeared frequencies of varieties harbored Hap2 (elite alleles G and T) with longer FL (> the mean of accessions with non-elite allele, 28.50 mm) and higher FS (> the mean of accessions with non-elite allele, 29.06 cN•tex-1) were 100 and 72.7%, respectively, which was higher than that of varieties harbored only on a single elite allele (G or T, 77.9 or 61.9%), suggesting a favorable haplotype for selecting varieties with superior FL and FS. These FMs could be valuable for the high-throughput selection of superior materials by providing genotypic information in cotton breeding programs.

A Genome Wide Association Study Revealed Key Single Nucleotide Polymorphisms/Genes Associated With Seed Germination in Gossypium hirsutum L.

Fast and uniform seed germination is essential to stabilize crop yields in agricultural production. It is important to understand the genetic basis of seed germination for improving the vigor of crop seeds. However, little is known about the genetic basis of seed vigor in cotton. In this study, we evaluated four seed germination-related traits of a core collection consisting of 419 cotton accessions, and performed a genome-wide association study (GWAS) to explore important loci associated with seed vigor using 3.66 million high-quality single nucleotide polymorphisms (SNPs). The results showed that four traits, including germination potential, germination rate, germination index, and vigor index, exhibited broad variations and high correlations. A total of 92 significantly associated SNPs located within or near 723 genes were identified for these traits, of which 13 SNPs could be detected in multiple traits. Among these candidate genes, 294 genes were expressed at seed germination stage. Further function validation of the two genes of higher expression showed that Gh_A11G0176 encoding Hsp70-Hsp90 organizing protein negatively regulated Arabidopsis seed germination, while Gh_A09G1509 encoding glutathione transferase played a positive role in regulating tobacco seed germination and seedling growth. Furthermore, Gh_A09G1509 might promote seed germination and seedling establishment through regulating glutathione metabolism in the imbibitional seeds. Our findings provide unprecedented information for deciphering the genetic basis of seed germination and performing molecular breeding to improve field emergence through genomic selection in cotton.

GhENODL6 Isoforms from the Phytocyanin Gene Family Regulated Verticillium Wilt Resistance in Cotton.

Verticillium wilt (VW), a fungal disease caused by Verticillium dahliae, currently devastates cotton fiber yield and quality seriously, yet few resistance germplasm resources have been discovered in Gossypium hirsutum. The cotton variety Nongda601 with suitable VW resistance and high yield was developed in our lab, which supplied elite resources for discovering resistant genes. Early nodulin-like protein (ENODL) is mainly related to nodule formation, and its role in regulating defense response has been seldom studied. Here, 41 conserved ENODLs in G. hirsutum were identified and characterized, which could divide into four subgroups. We found that GhENODL6 was upregulated under V. dahliae stress and hormonal signal and displayed higher transcript levels in resistant cottons than the susceptible. The GhENODL6 was proved to positively regulate VW resistance via overexpression and gene silencing experiments. Overexpression of GhENODL6 significantly enhanced the expressions of salicylic acid (SA) hormone-related transcription factors and pathogenicity-related (PR) protein genes, as well as hydrogen peroxide (H2O2) and SA contents, resulting in improved VW resistance in transgenic Arabidopsis. Correspondingly, in the GhENODL6 silenced cotton, the expression levels of both phenylalanine ammonia lyase (PAL) and 4-coumarate-CoA ligase (4CL) genes significantly decreased, leading to the reduced SA content mediating by the phenylalanine ammonia lyase pathway. Taken together, GhENODL6 played a crucial role in VW resistance by inducing SA signaling pathway and regulating the production of reactive oxygen species (ROS). These findings broaden our understanding of the biological roles of GhENODL and the molecular mechanisms underlying cotton disease resistance.

Blood compatibility evaluations of two-dimensional Ti3C2Txnanosheets.

Two-dimensional nanomaterial Ti3C2Txis a novel biomaterial used for medical apparatus. For its application, biosafety serves as a prerequisite for their usein vivo. So far, no research has systematically reported how Ti3C2Txinteracts with various components in the blood. In this work, we evaluated the hemocompatibility of Ti3C2Txnanosheets which we prepared by HF etching. Effects of the concentration and size of Ti3C2Txon the morphology and hemolysis rate of human red blood cells (RBCs), the structure and conformation of plasma proteins, the complement activation, as well asin vitroblood coagulation were studied. In general, Ti3C2Txtakes on good blood compatibility, but in the case of high concentration (>30 µg ml-1) and 'small size' (about 100 nm), it led to the rupture of RBCs membrane and a higher rate of hemolysis. Meanwhile, platelets and complement were inclined to be activated with the increased concentration, accompanying the changed configuration of plasma proteins dependent on concentration. Surprisingly, the presence of Ti3C2Txdid not significantly disrupt the coagulation.In vitrocell culture, the results prove that when the Ti3C2Txconcentration is as high as 60 µg ml-1and still has good biological safety. By establishing a fuzzy mathematical model, it was proved that the hemocompatibility of Ti3C2Txis more concentration-dependent than size-dependent, and the hemolysis rate is the most sensitive to the size and concentration of the Ti3C2Tx. These findings provide insight into the potential use of Ti3C2Txas biofriendly nanocontainers for biomaterialsin vivo.

High-quality genome assembly and resequencing of modern cotton cultivars provide resources for crop improvement.

Cotton produces natural fiber for the textile industry. The genetic effects of genomic structural variations underlying agronomic traits remain unclear. Here, we generate two high-quality genomes of Gossypium hirsutum cv. NDM8 and Gossypium barbadense acc. Pima90, and identify large-scale structural variations in the two species and 1,081 G. hirsutum accessions. The density of structural variations is higher in the D-subgenome than in the A-subgenome, indicating that the D-subgenome undergoes stronger selection during species formation and variety development. Many structural variations in genes and/or regulatory regions potentially influencing agronomic traits were discovered. Of 446 significantly associated structural variations, those for fiber quality and Verticillium wilt resistance are located mainly in the D-subgenome and those for yield mainly in the A-subgenome. Our research provides insight into the role of structural variations in genotype-to-phenotype relationships and their potential utility in crop improvement.

Cotton GhSSI2 isoforms from the stearoyl acyl carrier protein fatty acid desaturase family regulate Verticillium wilt resistance.

Lipids are major and essential constituents of plant cells and provide energy for various metabolic processes. However, the function of the lipid signal in defence against Verticillium dahliae, a hemibiotrophic pathogen, remains unknown. Here, we characterized 19 conserved stearoyl-ACP desaturase family proteins from upland cotton (Gossypium hirsutum). We further confirmed that GhSSI2 isoforms, including GhSSI2-A, GhSSI2-B, and GhSSI2-C located on chromosomes A10, D10, and A12, respectively, played a dominant role to the cotton 18:1 (oleic acid) pool. Suppressing the expression of GhSSI2s reduced the 18:1 level, which autoactivated the hypersensitive response (HR) and enhanced cotton Verticillium wilt and Fusarium wilt resistance. We found that low 18:1 levels induced phenylalanine ammonia-lyase-mediated salicylic acid (SA) accumulation and activated a SA-independent defence response in GhSSI2s-silenced cotton, whereas suppressing expression of GhSSI2s affected PDF1.2-dependent jasmonic acid (JA) perception but not the biosynthesis and signalling cascade of JA. Further investigation showed that structurally divergent resistance-related genes and nitric oxide (NO) signal were activated in GhSSI2s-silenced cotton. Taken together, these results indicate that SA-independent defence response, multiple resistance-related proteins, and elevated NO level play an important role in GhSSI2s-regulated Verticillium wilt resistance. These findings broaden our knowledge regarding the lipid signal in disease resistance and provide novel insights into the molecular mechanism of cotton fungal disease resistance.

A large-scale genomic association analysis identifies a fragment in Dt11 chromosome conferring cotton Verticillium wilt resistance.

Verticillium wilt (VW) is a destructive disease that results in great losses in cotton yield and quality. Identifying genetic variation that enhances crop disease resistance is a primary objective in plant breeding. Here we reported a GWAS of cotton VW resistance in a natural-variation population, challenged by different pathogenicity stains and different environments, and found 382 SNPs significantly associated with VW resistance. The associated signal repeatedly peaked in chromosome Dt11 (68 798 494-69 212 808) containing 13 core elite alleles undescribed previously. The core SNPs can make the disease reaction type from susceptible to tolerant or resistant in accessions with alternate genotype compared to reference genotype. Of the genes associated with the Dt11 signal, 25 genes differentially expressed upon Verticillium dahliae stress, with 21 genes verified in VW resistance via gene knockdown and/or overexpression experiments. We firstly discovered that a gene cluster of L-type lectin-domain containing receptor kinase (GhLecRKs-V.9) played an important role in VW resistance. These results proved that the associated Dt11 region was a major genetic locus responsible for VW resistance. The frequency of the core elite alleles (FEA) in modern varieties was significantly higher than the early/middle varieties (12.55% vs 4.29%), indicating that the FEA increased during artificial selection breeding. The current developmental resistant cultivars, JND23 and JND24, had fixed these core elite alleles during breeding without yield penalty. These findings unprecedentedly provided genomic variations and promising alleles for promoting cotton VW resistance improvement.

Tissue-specific expression of GhnsLTPs identified via GWAS sophisticatedly coordinates disease and insect resistance by regulating metabolic flux redirection in cotton.

Cotton (Gossypium hirsutum) is constantly attacked by pathogens and insects. The most efficient control strategy is to develop resistant varieties using broad-spectrum gene resources. Several resistance loci harboured by superior varieties have been identified through genome-wide association studies. However, the key genes and/or loci have not been functionally identified. In this study, we identified a locus significantly associated with Verticillium wilt (VW) resistance, and within a 145.5-kb linkage disequilibrium, two non-specific lipid transfer protein genes (named GhnsLTPsA10) were highly expressed under Verticillium pathogen stress. The expression of GhnsLTPsA10 significantly increased in roots upon Verticillium dahliae stress but significantly decreased in leaves under insect attack. Furthermore, GhnsLTPsA10 played antagonistic roles in positively regulating VW and Fusarium wilt resistance and negatively mediating aphid and bollworm resistance in transgenic Arabidopsis and silenced cotton. By combining transcriptomic, histological and physiological analyses, we determined that GhnsLTPsA10-mediated phenylpropanoid metabolism further affected the balance of the downstream metabolic flux of flavonoid and lignin biosynthesis. The divergent expression of GhnsLTPsA10 in roots and leaves coordinated resistance of cotton against fungal pathogens and insects via the redirection of metabolic flux. In addition, GhnsLTPsA10 contributed to reactive oxygen species accumulation. Therefore, in this study, we elucidated the novel function of GhnsLTP and the molecular association between disease resistance and insect resistance, balanced by GhnsLTPsA10. This broadens our knowledge of the biological function of GhnsLTPsA10 in crops and provides a useful locus for genetic improvement of cotton.

A stable QTL qSalt-A04-1 contributes to salt tolerance in the cotton seed germination stage.

KEY MESSAGE: A stable QTL qSalt-A04-1 for salt tolerance in the cotton seed germination stage, and two candidate genes, GhGASA1 and GhADC2, that play negative roles by modulating the GA and PA signalling pathways, respectively, were identified. The successful transition of a seed into a seedling is a prerequisite for plant propagation and crop yield. Germination is a vulnerable stage in a plant's life cycle that is strongly affected by environmental conditions, such as salinity. In this study, we identified a novel quantitative trait locus (QTL) qRGR-A04-1 associated with the relative germination rate (RGR) after salt stress treatment based on a high-density genetic map under phytotron and field conditions, with LOD values that ranged from 6.65 to 16.83 and 6.11-12.63% phenotypic variations in all five environmental tests. Two candidate genes with significantly differential expression between the two parents were finally identified through RNA-seq and qRT-PCR analyses. Further functional analyses showed that GhGASA1- and GhADC2-overexpression lines were more sensitive to salt stress than wild-type Arabidopsis based on the regulation of the transcript levels of gibberellic acid (GA)- and polyamine (PA)- related genes in GA and PA biosynthesis and the reduction in the accumulation of GA and PA, respectively, under salt stress. Virus-induced gene silencing analysis showed that TRV:GASA1 and TRV:ADC2 were more tolerant to salt stress than TRV:00 based on the increased expression of GA synthesis genes and decreased H2O2 content, respectively. Taken together, our results suggested that QTL qRGR-A04-1 and its two harboured genes, GhGASA1 and GhADC2, are promising candidates for salt tolerance improvement in cotton.

Evolution, expression and functional analysis of cultivated allotetraploid cotton DIR genes.

BACKGROUND: Dirigent (DIR) proteins mediate regioselectivity and stereoselectivity during lignan biosynthesis and are also involved in lignin, gossypol and pterocarpan biosynthesis. This gene family plays a vital role in enhancing stress resistance and in secondary cell-wall development, but systematical understanding is lacking in cotton. RESULTS: In this study, 107 GbDIRs and 107 GhDIRs were identified in Gossypium barbadense and Gossypium hirsutum, respectively. Most of these genes have a classical gene structure without intron and encode proteins containing a signal peptide. Phylogenetic analysis showed that cotton DIR genes were classified into four distinct subfamilies (a, b/d, e, and f). Of these groups, DIR-a and DIR-e were evolutionarily conserved, and segmental and tandem duplications contributed equally to their formation. In contrast, DIR-b/d mainly expanded by recent tandem duplications, accompanying with a number of gene clusters. With the rapid evolution, DIR-b/d-III was a Gossypium-specific clade involved in atropselective synthesis of gossypol. RNA-seq data highlighted GhDIRs in response to Verticillium dahliae infection and suggested that DIR gene family could confer Verticillium wilt resistance. We also identified candidate DIR genes related to fiber development in G. barbadense and G. hirsutum and revealed their differential expression. To further determine the involvement of DIR genes in fiber development, we overexpressed a fiber length-related gene GbDIR78 in Arabidopsis and validated its function in trichomes and hypocotyls. CONCLUSIONS: These findings contribute novel insights towards the evolution of DIR gene family and provide valuable information for further understanding the roles of DIR genes in cotton fiber development as well as in stress responses.

Dynamic characteristics and functional analysis provide new insights into long non-coding RNA responsive to Verticillium dahliae infection in Gossypium hirsutum.

BACKGROUND: Verticillium wilt is a widespread and destructive disease, which causes serious loss of cotton yield and quality. Long non-coding RNA (lncRNA) is involved in many biological processes, such as plant disease resistance response, through a variety of regulatory mechanisms, but their possible roles in cotton against Verticillium dahliae infection remain largely unclear. RESULTS: Here, we measured the transcriptome of resistant G. hirsutum following infection by V. dahliae and 4277 differentially expressed lncRNAs (delncRNAs) were identified. Localization and abundance analysis revealed that delncRNAs were biased distribution on chromosomes. We explored the dynamic characteristics of disease resistance related lncRNAs in chromosome distribution, induced expression profiles, biological function, and these lncRNAs were divided into three categories according to their induced expression profiles. For the delncRNAs, 687 cis-acting pairs and 14,600 trans-acting pairs of lncRNA-mRNA were identified, which indicated that trans-acting was the main way of Verticillium wilt resistance-associated lncRNAs regulating target mRNAs in cotton. Analyzing the regulation pattern of delncRNAs revealed that cis-acting and trans-acting lncRNAs had different ways to influence target genes. Gene Ontology (GO) enrichment analysis revealed that the regulatory function of delncRNAs participated significantly in stimulus response process, kinase activity and plasma membrane components. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis indicated that delncRNAs participated in some important disease resistance pathways, such as plant-pathogen interaction, alpha-linolenic acid metabolism and plant hormone signal transduction. Additionally, 21 delncRNAs and 10 target genes were identified as being involved in alpha-linolenic acid metabolism associated with the biosynthesis of jasmonic acid (JA). Subsequently, we found that GhlncLOX3 might regulate resistance to V. dahliae through modulating the expression of GhLOX3 implicated in JA biosynthesis. Further functional analysis showed that GhlncLOX3-silenced seedlings displayed a reduced resistance to V. dahliae, with down-regulated expression of GhLOX3 and decreased content of JA. CONCLUSION: This study shows the dynamic characteristics of delncRNAs in multiaspect, and suggests that GhlncLOX3-GhLOX3-JA network participates in response to V. dahliae invasion. Our results provide novel insights for genetic improvement of Verticillium wilt resistance in cotton using lncRNAs.

A high-density genetic map and multiple environmental tests reveal novel quantitative trait loci and candidate genes for fibre quality and yield in cotton.

KEY MESSAGE: A high-density linkage map of an intraspecific RIL population was constructed using 6187 bins to identify QTLs for fibre quality- and yield-related traits in upland cotton by whole-genome resequencing. Good fibre quality and high yield are important production goals in cotton (Gossypium hirsutum L.), which is a leading natural fibre crop worldwide. However, a greater understanding of the genetic variants underlying fibre quality- and yield-related traits is still required. In this study, a large-scale population including 588 F7 recombinant inbred lines, derived from an intraspecific cross between the upland cotton cv. Nongdamian13, which exhibits high quality, and Nongda601, which exhibits a high yield, was genotyped by using 232,946 polymorphic single-nucleotide polymorphisms obtained via a whole-genome resequencing strategy with 4.3-fold genome coverage. We constructed a high-density bin linkage map containing 6187 bin markers spanning 4478.98 cM with an average distance of 0.72 cM. We identified 58 individual quantitative trait loci (QTLs) and 25 QTL clusters harbouring 94 QTLs, and 119 previously undescribed QTLs controlling 13 fibre quality and yield traits across eight environments. Importantly, the QTL counts for fibre quality in the Dt subgenome were more than two times that in the At subgenome, and chromosome D02 harboured the greatest number of QTLs and clusters. Furthermore, we discovered 24 stable QTLs for fibre quality and 12 stable QTLs for yield traits. Four novel major stable QTLs related to fibre length, fibre strength and lint percentage, and seven previously unreported candidate genes with significantly differential expression between the two parents were identified and validated by RNA-seq. Our research provides valuable information for improving the fibre quality and yield in cotton breeding.